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CN Abstract 

The specific shear viscosity r) of a classically rotating system of nucleons that interact via a 

+-> 

monopole pairing interaction is calculated including the effects of thermal fluctuations and coupling 

O 

to pair vibrations within the selfconsistent quasiparticle random-phase approximation. It is found 
that J) increases with angular momentum M at a given temperature T. In medium and heavy 
J> systems, r\ decreases with increasing T at T > 2 MeV and this feature is not affected much by 

angular momentum. But in lighter systems (with the mass number A < 20), r\ increases with T at 
a value of M close to the maximal value M max , which is defined as the limiting angular momentum 
for each system. The values of n obtained within the schematic model as well as for systems 

(N 

y— I with realistic single-particle energies are always larger than the universal lower-bound conjecture 

h/(Ank B ) up to T=5 MeV. 
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I. INTRODUCTION 



Viscosity is an important mechanical property of fluids, which describes their resistance 
to flow. The recent ultrarelativistic Au-Au and Pb-Pb collisions at the Relativistic Heavy 
Ion Collider (RHIC) at Brookhaven National Laboratory pQ and the Large Hadron Collider 
(LHC) at CERN [2] have revealed a strongly interacting matter that behaves like a nearly 
perfect fluid with extremely low viscosity. This has generated a high interest in the study 
of viscosity in various systems, which consist of strongly interacting particles. In nuclear 
physics, although viscosity has been theoretically calculated and experimentally extracted 
in 1970s, the calculations of the specific shear viscosity fj = rj/s, that is the ratio of shear 
viscosity rj to the entropy volume density s, in finite nuclei as a function of temperature T 
was reported very recently only in two papers [21 H] • The results of these works, which have 
been carried out in two different approaches, show that the specific shear viscosity in hot 
nuclei at a temperature T as high as 5 MeV is actually very close to what obtained in the 
strongly interacting matter discovered at RHIC and LHC and T > 170 MeV. The approach 
in Ref. [I] used the Green-Kubo relation to calculate the shear viscosity rj of a finite hot 
nucleus directly from the width and energy of the giant dipole resonance (GDR) of this 
nucleus. The Green-Kubo relation expresses the shear viscosity in terms of the correlation 
function of the shear stress tensors. As a result, the shear viscosity can be calculated from 
the imaginary part of the retarded Green's function that describes the transport process. 

The studies in Refs. P, H] consider systems at zero angular momentum. How the shear 
viscosity changes in a hot rotating finite system is an interesting question. Recently the 
approach in Ref. [I] has been extended to calculate the specific shear viscosity from the 
parameters of the giant dipole resonance (GDR) in hot rotating nuclei |5j. Following the 
same line, in the present paper we would like to study the specific shear viscosity of a finite 
system of paired fermions, which interact via the monopole pairing force at finite angular 
momentum. We first consider a schematic model with doubly folded equidistant levels, then 
apply the formalism to several realistic nuclei. Despite being illustrated by the results of 
calculations within nuclei, the present formalism can be applied to any finite systems of 
fermions with discrete single-particle energies interacting via a monopole pairing force. 

The paper is organized as follows. The formalism is derived in Sec. [Hj The analysis of 
numerical results obtained in the schematic multilevel model with pairing as well as by using 
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realistic single-particle spectra for 20 O, 44 Ca, and 120 Sn are discussed in Sec. Ill The paper 
is summarized in the last section, where conclusions are drawn. 



II. FORMALISM 

We consider the pairing Hamiltonian describing a spherical system, which is classically 
rotating about the symmetry z axis [6]: 

H = Hp — XN — -yM . (1) 

Here A and 7 are the chemical potential and rotation frequency, respectively. Hp describes 
a system of N fermions, which interact via a monopole pairing force with the constant 
parameter G, namely 

Hp = y]e k (a\_ k a +k + a ] _ k a_ k ) - G^2a\a ] _ k a„ k ia kl , (k and k' > 0) , (2) 

k kk' 

where a} ±k (a± k ) are the creation (annihilation) operators of a nucleon (neutron or proton) 
with angular momentum k, projection ±m k (m k > 0), and energy e k . The particle-number 
operator N and total angular momentum M, which coincides with its z projection, are given 

as 

N = yXa^a+k + ai k a- k ) , M = ^2,m k {a\ k a +k - a]_ k a- k ) ■ (3) 

k k 

By using the Bogolyubov transformation a* k = u k a\ + v k a^ k from the particle operators, 
a k and a k to the quasiparticle ones, a k and a k , the Hamiltonian ([T]) is transformed into the 
quasiparticle one "H, whose explicit form can be found, e.g., in Refs. |7, 8]. 

The Hamiltonian ([T]) has been recently solved within the FTBCS1+SCQRPA in Refs. [3 
|S]. The FTBCS1+SCQRPA includes the effect caused by the quasiparticle-number fluc- 
tuation (QNF) at finite temperature (FT), which is neglected in the standard BCS so- 
lution of the pairing problem, as well as the effect due to coupling to pairing vibration 
within the selfconsistent quasiparticle random-phase approximation (SCQRPA). Since the 
FTBCS1+SCQRPA has already been discussed thoroughly in Refs. [1H9], we summarize 
below only their main results, which are necessary for calculations in the present paper. 
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A. FTBCS1 equations at finite angular momentum 



The FTBCS1 equations are a set of coupled equations for the pairing gap, particle number 
and angular momentum. The abbreviation "FTBCS1" denotes that, different from the con- 
ventional FTBCS, QNF is taken account in FTBCS1. The equation for the level-dependent 
pairing gap A k is given as 

A k = A + SA k , (4) 
with the level-independent gap A and the level-dependent gap 5A k defined as 

A = Gy^uyVh'^ - n>t' ~ n k') i £A fc = G + k _ _ u k v k , (5) 



1/ e k -Gvl-X\ 1/ _e k -Gvl-X 



E k / ' " 2 \ E 



where 



E k = ^(e k -Gvl-\y + Al, (6) 

and 5J\f k is the QNF at nonzero angular momentum 

5N% = (<W+) 2 + {5M k f = n+(l - n+) + n£ (1 - . (7) 

The equations for the particle number iV and total angular momentum M are [See Eqs. (25) 
and (26) of Ref. g]]: 

^ = 2]T vl{l-n+-n- k )+ l -{n + k +n- k ) , M = ]T m k {n+ - n~) . (8) 

k L -I & 

In these equations are quasiparticle occupation numbers at temperature T and angular 
momentum projection ^m k given by the Fermi-Dirac distribution of noninteracting fermions 

U±k = exp[(3(E kT im k )} + l ' = f ' (9) 
By neglecting the QNF, i.e., setting 5J\f k = in Eq. (|7|, the FTBCS1 equations become 
the conventional FTBCS ones. 

B. FTBCS1+SCQRPA 

Coupling to pair vibrations beyond the FTBCS1 is carried out by solving the SCQRPA 
equations for the pair vibration generated by the phonon operators 

el = £^=^=^. G^WH'- do) 

I v 1 - K - n k 
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The quasiparticle Hamiltonian H is then represented in the effective form as [8] 



k — k k' i-i kfj, 

where J\f k = c^ ±k a±k , A4 = M k + + A/jT, w M are the phonon energies (eigenvalues of the 
SCQRPA equations) and the vertex V k is given as 



(12) 



The explicit expressions of 6^, gk(k') and g^/ are given in Eq. (10) of Ref. [8], Eqs. (12) 
and (14) of Ref. [7j, respectively. The set of FTBCS1 and SCQRPA equations are solved 
selfconsistently to define the pairing gap Ajt, chemical potential A, rotation frequency 7, 
phonon energies u^, phonon amplitudes X k and y£ as well as the quasiparticle and phonon 
occupation numbers and at each values of temperature T and angular momentum M. 
After solving the FTBCS1+SCQRPA equations, the Green's function for the quasiparticle 



propagation is found by using the effective Hamiltonian (11) as 

G ±k (E) ' ' 



2ir E - E± - M±(E) ' 



(13) 



with the modified quasiparticle energies 



E k = h k + Ikk , 



(14) 



and the mass operator, which has the analytic continuation into the complex energy plane 
as 

M±(u ± is) = M±( U ) T . (15) 

The real and imaginary parts of this analytic continuation define, respectively, the energy 
shift and damping of the quasiparticle due to coupling to the SCQRPA, namely 



M\2 



[l-nt + iy^u-E 



(uj-E k 



+ £ 2 



1 - tit + 



n k +^)(u-E k + 



(uj-Et-Lorf + e* ' (u-Et + ^y + e* 



(16) 



(17) 



In Eqs. (15) - (17), v. 



QtQfi) is the phonon occupation number, which is found 



by solving a set of coupled equations (31) - (35) in Ref. [8], and e is a sufficient small 
parameter. 



The spectral intensities J k {oS) are derived from the relations [IT] 

J fc ± (o;)[exp(a;/r) + l]/2 = i[G ±k {u + ie) - G ±k {u - ie)\/2 = -lmG ±k (u) , (18) 
in the form 

and, the quasiparticle occupation numbers nf are found as 

/oo 
j£{u>)du . (20) 
-oo 

In the limit of small quasiparticle damping r y k (oj) — > 0, can be approximated with the 
Fermi-Dirac distribution 



+ 



1 



nl * -~— - , (21) 



exp(E fc ± /T) + 1 

where E^ are the solutions of the equations for the poles of the quasiparticle Green's func 



tions G± k (E) (13), namely 

Et-Et-Mf{Et) = 0. (22) 

A simplified version of FTBCS1+SCQRPA is called the FTBCS1+QRPA. In the latter, 
by ignoring the selfconsistency, the FTBCS1 is solved first. The obtained u k , v k , FTBCS1 
quasiparticle energies E k and rotation frequency 7 are then used to solve the QRPA equations 
to determine the phonon energies as well as the amplitudes and y^. The quasiparticle 



occupation numbers are approximated with the Fermi-Dirac distribution (21), where E k T 
"fink are used in place of E^ (See Eq. (24) of Ref. [8j), and the phonon occupation numbers 
are approximated with the Bose-Einstein distribution = [exp(o; At /T) — l] -1 . 



C. Specific shear viscosity 

As has been mentioned in the Introduction, the specific shear viscosity is defined as a 
function of T and M as 

V(T, M) = , s(T, M) = ^S(T, M) , (23) 

where 77 (T, M) is the shear viscosity and s(T, M) is the entropy density at finite T and M, 
p is the nuclear density, A is the number of nucleons, and S(T, M) is the entropy of the 
system. From Eq. (3) of Ref. [4], we have 

T){T,M) = 770— — — — — — . 24 

hm^o ImG(o;,0,0) 



A normalization is carried out in Eq. (24) so that in the limit of T = and M = 0, the value 



of r)(T, M) is equal to the parameter r] , whose value can be extracted from experimental 
systematics. 



By using Eq. (18), we obtain from Eq. (24) 



T)(T, M) = r)o 



J(0,T,M) 
J(0,0,0) 



(25) 



In Eq. (25) we introduce the shorthand notation 



J(0, T,M) = J2 4(0, T, M)/Q , J fc (0, T, M) = J+(0, T, M) + J" (0, T, M) , (26) 
where f2 is the sum of all single-particle levels]^] Replacing the right-hand side of Eq. (25) 



with its explicit expression obtained by using Eq. (19) with the mass operator (16) and 



quasiparticle damping (17), we obtain the final expression for rj(T, M) in the following form 

J+(T, M) + J-(T, M) 



" (T - M) ="»- 2/(0,0) 
with the explicit expressions for J^T, M) and J(0, 0) given as 



(27) 



J±(T,M) = J2 



+ 



l—n k 

(E±+o; Al )2+ e 2 1 (£±- ttV )2+ £ 2 



m 2 +eH em) 



J(0,0) = ^ 



(VP 2 



2 ' 



(28) 



6fc(0)+<? fcfe (0)+^(0) 



|} fc (0)+<? fcfc (0)+w M (0)] +e 



2 • 



(29) 



where 6^(0), ^(0), £"^(0), and w M (0) are, respectively, the values of bk, qtk-, solutions of Eq. 



(22) and SCQRPA eigenvalues, determined at T = and M — 0. 



The entropy S(T, M) is calculated as the sum of the quasiparticle and phonon entropies 



as 



(30) 



S(T,M) = S a (T,M) + S Q (T,M) , 
where the quasiparticle entropy S a (T,M) is given as 

S a (T, M) = -J2 H In nt + (1 - n+) In(l - n+) + Inn" + (1 - n") ln(l - n")] , (31) 



1 For the spherical case: J(0,T,M) = V ,<>,./,;(). 7 . .W i V () ; with % =j+ 1/2. 
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whereas the phonon entropy is defined as the boson one, namely 

S Q (T,M) = ^[(l + ^)ln(l + ^)-^lni/J . (32) 

Within the FTBCS (FTBCS1), the mass operator M k (u), quasiparticle damping 7 & (o>), 
phonon occupation numbers and phonon energies are zero. Therefore, the expressions 
for J ± (T, M) and J(0, 0) are given as 

where E k = E^^L^m^ with E^. being the quasiparticle energies obtained within the FTBCS 
(FTBCS1). 

III. ANALYSIS OF NUMERICAL RESULTS 

The numerical calculations are carried out for the schematic model as well as realistic 
nuclei. The schematic model consists of Q doubly folded equidistant levels, which interact 
via a pairing force with parameter G. Before switching on the pairing interaction, Q/2 
lowest levels are occupied by N = Q particles with 2 particles on each levels with spin 
projections im/c = ±1/2, ±3/2, ±fi — 1/2, whereas Q/2 upper levels are empty. The 
distance between levels is chosen equal to 1 MeV [8]. As for the realistic nuclei, 20 O, 
44 Ca, and 120 Sn are considered. The single-particle spectra for these nuclei are obtained 
within the axially deformed Woods-Saxon potential including the spin-orbit and Coulomb 
interactions |10j . All the bound (negative energy) single-particle states are taken into account 
in the calculations. The pairing interaction parameter G is adjusted so that the pairing gap 
at T = reproduces the experimental value obtained from the odd-even mass difference. 
These nuclei have the proton closed shell, so there is only neutron pairing gap. Therefore 
G N is chosen equal to 1.04 MeV, 0.53 MeV, and 0.14 MeV, which yields Ajv(0) « 3.0 MeV, 
2.0 MeV, and 1.42 MeV for 20 O, 44 Ca, and 120 Sn, respectively. Since the pairing gap Q 
is level-dependent, in the analysis we plot the level- weighted gap A, which is defined as 
A = Afc/f2 with Q being the total number of levels. 

In the calculations of the shear viscosity rj we adopt the value of nuclear density p = 



0.16 fm . A value e = 0.5 MeV in Eqs. (16) and (17) is used in all calculations to ensure 



the smooth behavior of the obtained r/(T, M). Regarding the value 77(0, 0), in principle, a 



superfluid has zero viscosity at T = and M = 0. At T ^ (or M ^ 0), the pairing gap 
decreases with increasing T (or M) whereas the admixture of normal mode appears because 
of the QNF. As a result the shear viscosity r)(T, M) becomes finite at T ^ (or 0). As 
for the realistic nuclei under consideration, they consist of proton closed shells and neutron 
opened shells. Therefore the shear viscosity is finite even at T = because of the protons 
in the normal phase. Meanwhile, according to the lower bound f]/s > fiji^nkB) conjectured 
in Ref. [12] . the shear viscosity rj(T, M) cannot be at T ^ even at very small T because 
s is always positive at any T ^ 0. In other words, there are no perfect fluids, that is fluids 
having no shear viscosity, at T ^ 0. Since in the present paper we are interested in the 
behavior of 1](T, M) as a function of T and M, we adopt rj(0, 0) = r)o — 1.0 x 10 -23 Mev s 
fm~ 3 in the calculations for realistic nuclei as well as for the schematic model. This value 
is obtained by fitting the width of giant resonances at T = (See Ref. [I] for the detail 
discussions on the parameter r] ). 

Shown in Figs. [l]-[3]are the level- weight pairing gap A, density entropy s(T,M), shear 
viscosity 77 (T, M), and specific shear viscosity fj(T, M) = r](T, M) /s(T, M) obtained as func- 
tions of T at 3 values of angular momentum M for the schematic model with N = 10, 20 and 
100, respectively. As has been discussed in Refs. [8] and references therein, a distinguished 
feature of the FTBCS1 is that the gap A never collapses as the FTBCS gap at the critical 
temperature T c of the phase transition from the superfluid phase to the normal one, but, be- 
cause of QNF, it decreases monotonically with increasing T, and has a tail up to T = 5 MeV 
[Panels (a) - (c)]. With increasing N, the tail becomes depleted and eventually vanishes 
at iV — > 00. The gap decreases with increasing M as well known within the conventional 
FTBCS. The corrections due to coupling to SCQRPA slightly increases the pairing gap in 
the temperature region around the BCS phase transition, but this effect quickly decreases 
with increasing N to become negligible at N > 100. 

The entropy densities predicted by the FTBCS, FTBCS1, and FTBCS1+SCQRPA are 
close to each other, although the SCQRPA corrections slightly decrease the entropy at T < 
3 MeV for N = 10, and increase it at all T for N > 20 [Panels (d) - (f)]. Meanwhile, the 
FTBCS1+QRPA entropy density differs strongly with these results at high T because the 
Bose-Einstein distribution, which approximates the phonon occupation number within 
the QRPA, differs significantly from obtained by solving the selfconsistent equations (31) 
- (35) in Ref. [8] employed within the SCQRPA (See e.g. Fig. (5) of Ref. (7]). By comparing 
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FIG. 1. (Color online) Level- weight pairing gap A, density entropy s(T,M), shear viscosity 
rj(T, M), and specific shear viscosity rj(T, M) = rj(T, M)/s(T, M) obtained for the schematic model 
with N = 10 within the FTBCS (dotted line), FTBCS1 (dashed line), FTBCS1+QRPA (dash- 
dotted line), and FTBCS1+SCQRPA (solid line) at different values of angular momentum M. The 
horizontal dashed lines in (j) - (1) denotes the lower bound conjecture H/(AirkB) proposed in Ref. 
H2l. 



the panels (d) - (f) in these figures, one can see that angular momentum has negligible effect 
on the entropy density. 

Within the FTBCS, the shear viscosity rj increases with T from its value 770 at T = to 
reach at T > T c a value larger than r/ by one to two orders of magnitude, demonstrating 
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FIG. 2. (Color online) The same as in Fig. [I] but for N=20. 



that the normal fluid is highly viscous and its viscosity remains essentially constant in 
the normal phase [Panels (g) - (i)]. QNF within the FTBCS1 significantly decreases rj at 
T > T c , whereas the effects due coupling to the QRPA and SCQRPA increase and decrease 
rj for N < 20 and N > 20, respectively. The largest difference between the FTBCS1 and 
FTBCS1+SCQRPA is almost 2 orders of magnitude at T ~ 2.7 MeV. Increasing the angular 
momentum leads to an overall increase of rj, especially at low T. In other words, rotation 
makes a systems of fermions that interact via monopole pairing force becomes more viscous. 
As the particle number N increases, rj obtained within the FTBCS1+SCQRPA (QRPA) 
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FIG. 3. (Color online) The same as in Fig. [I] but for N=100. 



decreases to approach the prediction by the FTBCS1 [See Fig. [2](g) -l2j(i)], which increases 
with N. In large systems, e.g. N =100, rj obtained within the FTBCS1+SCQRPA (QRPA) 
becomes even smaller than that predicted by the FTBCS1 [See Fig. |3). The combine 
effect of temperature and angular momentum leads to a local minimum in the temperature 
dependence of 77 at T slightly below T c and local maximum at T slightly above T c . With 
increasing N this feature transforms into a local maximum at T ~ T c with a slight local 
minimum at each side. 

The overall trend of the specific shear viscosity 7/ is a decrease with increasing T in the 
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FIG. 4. (Color online) The same as in Fig. 1 but for u O 
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region T < T c . At T around T c a local maximum is seen, which is the direct consequence 
of the local maximum in the temperature dependence of rj discussed above. At T ^> T c the 
specific shear viscosity r] obtained with the FTBCS1+SCQRPA for large N (> 100) also 
decreases with increasing T. 

The qualitative features obtained for the shear viscosity with the schematic model as 
a function of T at various M also hold for realistic nuclei, namely 20 O, 44 Ca, and 120 Sn, 
as shown in Figs. [4]-[6j In a light system such as 20 O, QNF and coupling to SCQRPA 
(QRPA) cause dramatic effects, which significantly reduce rj (T, M) as compared to the 
prediction within the FTBCS (FTBCS1) in the temperature region T c < T < 3.5 MeV. For 
44 Ca and 120 Sn, the shear viscosity rj always decreases with increasing T at all M except for 
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FIG. 5. (Color online) The same as in Fig. 1 but for 44 Ca. 



the prediction by the FTBCS1+QRPA at M = 2h for 44 Ca because of the boson occupation 
numbers discussed above. These figures show that, at T = 5 MeV and M =0, the specific 
shear viscosity fj in spherical systems of paired fermions with realistic single-particle energies 
ranges from between twice to ten times larger than its value at T = and M = 0. Classical 
rotation strongly increases the specific shear viscosity in the light system such as 20 O at hight 
T, whereas in the heavier systems such as 44 Ca and 120 Sn such increase is much weaker. The 
change in the temperature dependence of r] in 20 O at high T from decreasing as T increases 
at low M to increasing with T at large M shows that rotation of a light system at high 
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FIG. 6. (Color online) The same as in Fig. 1 but for 120 Sn. 



T drives it from the behavior of a liquid to that of a gas. Meanwhile, medium and heavy 
systems behave as liquids because the predictions by the FTBCS1+SCQRPA show that 
their shear viscosity always decreases as T increases at high T and any M. In general, it 
is clear to see that the specific shear viscosity obtained within the schematic model as well 
as realistic nuclei considered here is always higher than the lower bound conjecture (KSS 
limit) proposed in Ref. [12] up to T=5 MeV. 

For small N, FTBCS1+SCQRPA (QRPA) predicts an increase fj with T at the largest 
value of M, namely 8 and 12 h for N = 10 and 20, respectively. This largest value should 
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FIG. 7. (Color online) Shear viscosity rj(T, M) and quasiparticle energies E£(T, M) (l22b for several 
levels within the schematic model with N = 10 [(a) - (d)] and for neutrons in Ca [(d) - (f)]. 



be smaller than the value M max , at which at least one solution of Eq. (22) becomes at 



a certain T, and turns negative at higher T. This is demonstrated in Fig. [7| where 77 and 
E k for several levels within the schematic model for N = 10 and for neutrons in 44 Ca are 
plotted against T at three values of M. For N = 10, at M max = 10ft one can see that rj 
has a singularity at T ~ 4.7 MeV [Fig. [7] (a)], where E£ crosses [Fig. [7] (c)]. For 44 Ca, 
a similar feature takes place at T = 3 MeV and M max = 3h [Figs. [7] (d) and [7| (f)]. The 



quasiparticle energies E k remain positive at all T and M because E k in Eq. (14) are always 
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positive. The limiting values M max of the angular momentum, determined for N = 20, 100, 
20 O, and 120 Sn, are 13h, 43h, 3h and 4h, respectively. That's is the reason why only the 
results obtained for M < M max are shown in Figs. [Ifj6j 

It is important to keep in mind that the results and analyses in the present paper are ob- 
tained by using the Hamiltonian ([I]), which describes essentially the motion of quasiparticles 
in the quasiparticle mean field, modified by rotation and coupling to vibration of monopole 
quasiparticle pairs. In general, they are not valid for nuclear systems whose residual interac- 
tions consist of nonzero multipolarities such as dipole, quadrupole, octupole, etc, since the 



latter correspond to the Green's functions different from Eq. (13), which lead to different 
behaviors of the shear viscosity r](T,M). As a matter of fact, the specific shear viscosity 
fj extracted in Ref. jl] based on the parameters of GDR in hot rotating nuclei shows a 
tendency of decreasing with increasing M at high T. However, the present work provides 
a theoretical formalism, which is based on the spectral intensities or response functions to 
derive and calculate the shear viscosity for not only nuclear but also any finite hot rotating 
systems. 



IV. CONCLUSIONS 

In the present paper, the Green-Kubo relation is used to calculate the specific shear 
viscosity from the retarded Green's function that describes the propagation of quasiparticles 
within the quasiparticle mean field of a classically rotating system of nucleons that interact 
via a monopole interaction. Thermal fluctuations such as the QNF are included within the 
FTBCS1, whereas coupling to monopole pair vibrations is taken into account within the 
SCQRPA. 

The general feature of the specific viscosity fj of this system can be summarized as follows. 
At a given temperature T, fj increases with the angular momentum M, that is a rotating 
system of paired fermions is more viscous. In medium and heavy systems, fj decreases with 
increasing T at T > 2 MeV and this feature is not affected much by angular momentum. 
However, in light systems, it increases with T at the values of angular momentum M close 



to M max , at which the quasiparticle energy E£ (22) crosses at a certain T. At T < 2 
MeV, local minima and/or local maximum appear because of the significant change in the 
curvature of the temperature dependence of the thermal pairing gap. Thermal fluctuations 
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and coupling to the quasiparticle pair vibrations within the SCQRPA significantly increases 
fj for small N systems with N < 10, whereas fj decreases for large iV > 10 systems. All the 
results of fj obtained within the schematic model as well as realistic nuclei are always larger 
than the universal lower bound of the specific shear viscosity up to T=5 MeV. 

Although the present paper considers the pairing model only, it lays the way of calculating 
the shear viscosity in not only nuclear but also any finite hot rotating Fermi systems based 
on the spectral intensities or response functions, and therefore opens the possibility to study 
the evolution of the shear viscosity as a function of temperature and angular momentum in 
these systems. 
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